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Abstract 



We present a study of inclusive direct-photon and pion production in hadronic 
interactions, focusing on a comparison of the ratio of 7/71"° yields with ex- 
pectations from next-to-leading order perturbative QCD (NLO pQCD). We 
also examine the impact of a phenomenological model involving kx smear- 
ing (which approximates effects of additional soft-gluon emission) on absolute 
predictions for photon and pion production and their ratio. 

PACS numbers: 12.38. Qk, 13.85.Qk,13.85.Ni 
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I. INTRODUCTION 



Direct-photon production in hadronic collisions at high transverse momenta (pr) has long 
been viewed as an ideal testing ground for the formalism of pQCD. A reliable theoretical 
description of the direct-photon process is of special importance because of its sensitivity to 
the gluon distribution in a proton through the quark-gluon scattering subprocess (gq — > jq). 
The gluon distribution, G(x), is relatively well constrained for x < 0.1 by data on deep- 
inelastic lepton-nucleon scattering (DIS) and dilepton ( "Drell-Yan" ) production, and by 
collider data on jet production at moderate x (0.1 — 0.25) [Q, but is less constrained at 
larger x 0. In principle, fixed-target direct-photon production can constrain G(x) at large x, 
and such data have therefore been incorporated in several modern global parton distribution 
function (PDF) analyses [0,g,(§. 

However, both the completeness of the theoretical NLO description of the direct-photon 
process, as well as the consistency of results from different experiments, have been subjects 
of intense debate fl]-|l0"fl. The inclusive production of hadrons provides a further means of 
testing the predictions of the NLO pQCD formalism. Deviations have been observed be- 
tween measured inclusive direct-photon and pion cross sections and NLO pQCD calculations. 
Examples of such discrepancies are shown in Figs. [I] and |2|, where ratios of data to theory 
are displayed as a function of xt = 2p T /y/s for p T > 3 GeV/c for photon and pion data, 
respectively (Unless otherwise indicated, all NLO calculations p3j-p9| in this paper 

use a single scale of fi = Pt/2, CTEQ4M PDFs and BKK fragmentation functions for 
pions [Q.) The uncertainties shown in the figures reflect mainly the statistical fluctuations 
in the measurements, except for some pioneering experiments from the ISR, which reported 
statistical errors combined with some of their systematic uncertainties. 

Several of the older experiments quote large (of order 40%) uncertainties in their normal- 
izations, making it difficult to judge the significance of the discrepancies observed in Figs. [I] 
and |2|. Nevertheless, it is clear that both the photon and pion data show disagreement 
with theory, and the different data sets also disagree. This is not surprising, because, espe- 
cially for direct-photon production, signals are often very small and difficult to extract, and 
backgrounds are quite large, especially at lower energies. Several experiments show better 
agreement with NLO pQCD than others, but the results do not provide great confidence in 
the theory nor in the quality of all the data. While it has been suggested that deviations 
from theory for both photons and pions can be ascribed to higher-order effects of initial-state 
soft-gluon radiation |],f7|, it seems unlikely that theoretical developments alone will be able 
to accommodate the level of scatter observed in Figs. [L] and 

II. 7/vr° RATIO 

Given the scatter of the data shown in Figs. [TJ and it may be instructive to consider 
measurements of the 7/71"° ratio in different experiments over a wide range of y/s. Both 
experimental and theoretical uncertainties tend to cancel in such a ratio, and this ratio 
should also be less sensitive to incomplete treatment of gluon radiation. In addition, since 
7T° production constitutes the primary background to photon production, the 7/71" ratio 
can serve as a measure of the difficulty of the experimental environment for making the 
direct-photon measurement. 
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The ratio of direct-photon to 7r° cross sections for both data and NLO calculations 
(solid curves) are shown, as a function of xt, in Figs. 0-0. Calculations for experiments 
using nuclear targets (E629, NA3, NA24, E706) were corrected approximately for nuclear 
dependence effects; these corrections were < 20% for 7r°'s and < 10% for photons. The 
results from E704 0, E629 @, and NA3 0, all at V« = 19.4 GeV, are displayed in 
Fig. |3|. For all three measurements, theory is high compared to data, somewhat less so 
for E629. Figure | shows the 7 to vr° ratio for WA70 0, NA24 |T§, and UA6 JTJ at 
y/s w 23 — 24 GeV. Just as for the comparisons at y/s = 19.4 GeV, theory is high relative to 
data, with the deviation between theory and data for WA70 and NA24 being greater than 
that for UA6. 

For larger values of y/s, the theoretically calculated value of the 7/71- ratio agrees better 
with experiment. Figure || shows the data from R806 (18j at y/s = 31 GeV, and E706 [|17j at 
y/s = 31.6 and 38.8 GeV. Theory provides a reasonable description of these measurements. 



For the results shown in Fig. o, theory is now somewhat below the data for R806 [18] at 



y/s = 45 and 53 GeV, but agrees with the R807 data |JT| at y/s = 53 GeV. Theory is also 
slightly low relative to data at y/s = 63 GeV for R806 R807 |20[, and R110 J21 



as shown in Fig. [7[ A compilation of these results, shown for simplicity without their 
uncertainties, is presented in Fig. || In this figure, the ratios of data to theory for the 7 to 
7T° measurements shown in Figs. [3|-[^] have been fitted to a constant value at high-pr, and the 
results plotted as a function of y/s. (Because the ratio depends somewhat on p T , we excluded 
from the fit points at low px, until a \ 2 /d.o.f. = 1.5 was reached for each fit.) The figure 
suggests an energy dependence in the ratio of data to theory for 7/71- production. There is, 
however, an indication of substantial differences between the experiments at low y/s, where 
the observed 7/71° is smallest, which makes it difficult to quantify this trend. Recognizing 
the presence of these differences is especially important because thus far only the low energy 
photon experiments have been used in PDF fits to extract the gluon distribution. 



III. PHENOMENOLOGICAL CORRECTIONS FOR SOFT GLUON EMISSION 

In hadronic hard-scattering processes, there is generally a substantial amount of effective 
parton transverse momentum, kj>, in the initial state resulting from gluon emission |7|,|32 



The presence of k^ impacts the final state and has been observed in measurements of Drell- 
Yan, diphoton, and heavy quark production; the amount of k^ expected from NLO calcu- 
lations is not sufficient to describe the data. The effective values of (kx) /parton for these 
processes vary from w 1 GeV/c at fixed target energies to 3-4 GeV/c at the Tevatron Col- 
lider. The growth is approximately logarithmic with center-of-mass energy J7,3"2"|. The size 



of the (kr) values, and their dependence on energy, argue against a purely "intrinsic" non- 
perturbative origin. Rather, the major part of this effect is generally attributed to soft 
gluon emission. While the importance of including gluon emission through the resummation 
formalism has long been recognized and such calculations have been available for some time 
for Drell-Yan fl3H|1 , diphoton p4| , |35| , and W/Z production |54| , they have only recently been 



developed for inclusive direct-photon production |36|-[4T|. 

In the absence of a more rigorous theoretical treatment of the impact of soft-gluon 
emission on high-p<r inclusive production, a more intuitive, but successful, phenomenological 
approach has been utilized for describing the effects of soft-gluon radiation f7[. The soft- 
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gluon radiation was parametrized in terms of an effective (kx) that provided an additional 
transverse impulse to the outgoing partons. Because of the steeply falling cross section in 
Px, such a (kx) can shift the production of final-state particles from lower to higher values 
of pt, effectively enhancing the cross section. Using this intuitive picture, the effects of soft 
gluon emission are approximated through a convolution of the NLO cross section with a 
Gaussian kx smearing function. The value of (kx) used for each kinematic regime is based 
on experimental data M. 

As described in Ref. [0], a leading order (LO) pQCD calculation |42| is used to generate 
correction factors (ratios of calculations for any given (kx) to the result for (kx) = 0) for 
inclusive cross sections (Fig. |). These pr-dependent K-factors are then applied to the NLO 
pQCD calculations. This procedure involves a risk of double-counting since some of the 
/cT-enhancements may already be contained in the NLO calculation. However, the effects of 
such double-counting are expected to be small @ . The relationship between this phenomeno- 
logical fc^-smearing and the Collins-Soper-Sterman (CSS) resummation formalism [13] has 



been discussed in Ref. 44 



As illustrated in the upper part of Fig. |9|, the K-factors for direct-photon production in 
E706 are large over the full range of pt, and have pr-dependent shapes. The lower part of 
Fig. |S] displays K-factors for 7r° production, based on the same model. For the same values 
of (kx), the K-factors in ir° production are somewhat smaller than in 7 production. This 
is reasonable because 7r° mesons originate from the fragmentation of partons, and therefore 
carry only a fraction of the kx given to the jet. Because the momentum fraction of the 
jet carried by a 7r° tends to increase with increasing px of the 7r° (sometimes referred to as 
trigger-bias), the kx corrections in the case of 7r° production are smaller and tend to have a 
weaker p^-dependence than the corrections for direct photons. 

In this model, the size of the correction depends sensitively on the value used for (kx)- 
Changes of 200 MeV/c produce substantial differences in the size of the correction. Values of 
(kx) are only accurate to about this range J7J, and so, even within the context of this model, 
it is difficult to obtain the kind of precision needed for extracting global parton distributions. 
In addition, there are different ways to implement the fc^-corrections |4[] , which can produce 
quantitative differences in the magnitude and shape of the kx correction factors. 

We have investigated the kind of K-factors that would be expected for direct-photon 
production from parton-showering models |45];|46| . These programs do not provide sufficient 
smearing at fixed-target energies because shower development is constrained by cut-off pa- 
rameters that ensure the perturbative nature of the process. Consequently, these calculations 
allow additional input kx for Gaussian smearing, and are often used that way in comparisons 
to data (default values are (k T ) = for herwig and (k T ) = 0.9 GeV/c for pythia). Using 
default settings for other program parameters and an input (kx) of 1.2 GeV/c (to compare 
to our previous results 0) for the smearing, relative to these same settings with kx = in 



these generators, we obtain the K-factors shown in Fig. 10. These corrections are different 



than those found using the LO pQCD prescription [42 . These differences must be kept in 



mind when making comparisons to data. Applying the K-factors in Fig. TD to NLO pQCD 
results in the comparison to data from E706 shown in Fig. |TI[ 

The dashed curves in Figs. describe the predicted ratios of 7 and n° cross sections 
using the previous kx corrections JFj. The impact of kx corrections upon 7/7?° is generally 
minimal, indicating that the trend observed in Fig. IS] cannot be easily understood on the 
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basis of corrections for kx alone. 



IV. RECENT PROGRESS IN THEORY 



Resummed pQCD calculations for single direct-photon production are currently under 
development f36H4T|. Substantial corrections to fixed-order QCD calculations are expected 



from soft-gluon emission, especially in regions of phase space where gluon emission is re- 
stricted kinematically. Such restrictions occur, for example, in the production of Drell-Yan 
or diphoton pairs at low p?, where effects of soft-gluon emission have to be resummed in 
order to obtain adequate description of cross sections f33| . In addition, at high x, there is 
a similar suppression of gluon radiation due to the rapidly falling parton distributions. A 
complete description of the cross section in this region requires the resummation of "thresh- 
old" terms. For inclusive "single-arm" observables, there is no restriction in phase space 
for soft gluon radiation, and therefore no enhancement to the cross section is expected 
from traditional resummation calculations except at high x |[47|| . Two recent independent 



threshold-resummed pQCD calculations for direct photons exhibit far less depen- 

dence on QCD scales than found in NLO theory. These calculations agree with the NLO 
prediction for the scale \x ~ Pt/2 at low pt (without inclusion of explicit hr or recoil effects) , 
and show an enhancement in cross section at high p? due to the threshold terms. 

A method for simultaneous treatment of recoil and threshold corrections in inclusive 
single-photon cross sections is being developed within the formalism of collinear fac- 



torization. This approach accounts explicitly for the recoil from soft radiation in the 
hard-scattering subprocess, and conserves both energy and transverse momentum for the 
resummed radiation. The possibility of substantial enhancements from higher-order per- 
turbative and power-law nonperturbative corrections relative to NLO are indicated at both 
moderate and high pt for fixed-target energies, similar to the enhancements obtained with 
the simple fc^-smearing model discussed above. 

Figure [12| displays the results of an example calculation based on the latest approach 
involving both threshold and recoil resummation for comparison with direct-photon mea- 
surements from E706 [pfj . As shown in Fig. O, this theoretical result is substantially higher 



than the prediction from NLO pQCD, higher than the theory using just threshold resum- 



mation, and closer to the phenomenological A^-smearing model used in Fig. 11 



V. CONCLUSIONS 

We have examined experimental information on the production of direct photons and 
pions at large px and reviewed how LO pQCD has been used to estimate the impact of 
k? on their inclusive production. Simple phenomenological models can improve agreement 
between pQCD calculations and much of the inclusive data over a wide range of ^Js [[?]. 
While the deficiencies of such models are clear, and have been discussed in recent literature, 
the emerging formalism for the full (threshold and recoil) resummation of inclusive direct- 



photon cross sections [41] appears to support some of the understanding of soft multiple 



gluon emission that has been gained using the approximate kr tools. 
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Despite the apparent inconsistencies between different direct-photon experiments |7| -||. 
we found it instructive to consider results on the measured 7/71" ratios since various exper- 
imental and theoretical uncertainties tend to cancel in such ratios. We find that the ratios 
taken from theory agree reasonably with those computed from data for y/s ^ 30 GeV. How- 
ever, there appears to be a systematic trend with energy that is unexplained, and differences 
between experiments that are particularly signficant at lower energy. 

While there is still no resummation calculation for inclusive pion production, the trend 
of recent developments in direct-photon processes has led to an increased appreciation of 
the importance of the effect of multiple gluon emission, and to the emergence of tools for 
incorporating these effects. These latest theoretical developments encourage optimism that 
the long-standing difficulties in developing an adequate description of these processes can 
eventually be resolved, making possible a global re-examination of parton distributions with 



an emphasis on the determination of the gluon distribution from the direct-photon data 48 
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FIG. 1. Comparison between proton-induced direct-photon data and NLO pQCD calculations 
as a function of photon xt- 
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FIG. 2. Comparison between proton-induced pion data and NLO pQCD calculations for ex- 
periments spanning yfs = 19.4 to 900 GeV as a function of xt- 
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FIG. 3. Comparison of 7/V rates as a function of xt for E704, E629, and NA3 at 
y/s = 19.4 GeV. Overlayed are the results from NLO pQCD (solid) and fc^-enhanced calcula- 
tions (dashed). Values of (kx) used for the fey-enhanced calculations are given in the legend. 
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FIG. 4. Comparison of 7/71" rates as a function of xt for WA70, NA24, and UA6 at 
yfs ~ 23 — 24 GeV. Overlayed are the results from NLO pQCD (solid) and /cr-enhanced cal- 
culations (dashed). Values of (hp) used for the fc^-enhanced calculations are given in the legend. 
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FIG. 5. Comparison of 7/vr rates as a function of xt for R806 at \fs = 31 GeV and E706 
at y/s = 31.6 and 38.8 GeV. Overlayed are the results from NLO pQCD (solid) and /cr-enhanced 
calculations (dashed). Values of (kx) used for the fey-enhanced calculations are given in the legend. 
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FIG. 6. Comparison of 7/vr rates as a function of xt for R806 at \fs = 45 and 53 GeV 
and R807 at yfs = 53 GeV. Overlayed are the results from NLO pQCD (solid) and /cr-enhanced 
calculations (dashed). Values of (kx) used for the /^-enhanced calculations are given in the legend. 
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FIG. 7. Comparison of 7/71" rates as a function of for R806, R807, and R110 at 
= 63 GeV. Overlayed are the results from NLO pQCD (solid) and fe^-enhanced calculations 
(dashed). Values of (kx) used for the fey-enhanced calculations are given in the legend. 
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the direct-photon experiments considered in this paper. The values represent fits to the ratio of 
data to NLO pQCD theory, without /cr-enhancement (see text). 
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FIG. 9. The dependence of the fc^-enhancement on p?, K(pt), for variations in (kx) for the 
E706 pBe data at ^fs = 31.6 GeV. 
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FIG. 10. Comparison between the kx enhancements, K(pr), from the LO pQCD calculation 
(dashed), herwig (dotted), and pythia (dashed-dotted). 
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FIG. 11. Comparison between the E706 direct-photon data at y/s = 31.6 GeV and the NLO 
pQCD calculation (solid), and the NLO theory enhanced by K-factors associated with the LO 
calculation (dashed), herwig (dotted), and pythia (dashed-dotted). 
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FIG. 12. Direct-photon cross section for the E706 data at \fs = 31.6 GeV compared to recent 
QCD calculations. The dotted line represents the full NLO calculation [28], while the dashed and 
solid lines, respectively, incorporate purely threshold resummation [36| and joint threshold and 
recoil resummation ||41]|. 
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